o 



o 






X 



EPJ manuscript No. 

(will be inserted by the editor) 



Tendencies toward nematic order in YBa2Cu306+(^: 
Uniform distortion vs. incipient charge stripes 

Matthias Vojta 

Institut fiir Theoretische Physik, Universitat zu Koln, Ziilpicher Strafie 77, 50937 Koln, Germany 



Abstract. Recent neutron scattering and transport data obtained on underdoped 
^vj ' YBa2Cu306+5, with strong signatures of rotation symmetry breaking at low tem- 

^SJ , peratures, point toward electron-nematic order in the charge sector. Such order 

may originate from a uniform distortion with d-wave symmetry or as a precursor 
of a uni- directional stripe phase. Here, we discuss whether the neutron scattering 
data can be linked to incipient charge stripes. We employ and extend a phe- 
nomenological model for collective spin and charge fluctuations and analyze the 
resulting spin excitation spectrum under the influence of lattice anisotropies. Our 
results show that the experimentally observed temperature-dependent magnetic 
incommensurability is compatible with a scenario of incipient stripes, the temper- 



§^' 

^ ' ature dependence being due to the temperature variation of both strength and 

correlation length of the charge stripes. Finally, we propose further experiments 
C^ ■ to distinguish the possible theoretical scenarios. 



^ • 1 Introduction 

o 

^ O^ ■ The interplay of different ordering tendencies is a recurring theme in unconventional supercon- 

ductors - this applies to cuprates, iron pnictides, and heavy-fermion superconductors. In high-Tc 
superconducting cuprates, phenomena of spontaneous lattice symmetry breaking have been ob- 
^ ' served inside the pseudogap regime [I1I2]. In particular, stripe order in both the spin and charge 

^O ! sectors has been established in the La2-£cSr2;Cu04 (or 214) family 3 , while in YBa2Cu306+5 

^«0 ■ (YBCO), signatures of electronic nematic order have been reported |4il5j. Given the similarities 

in the spin excitation spectra at intermediate energies - the so-called hour-glass spectrum - 
across different cuprate families, it has long been proposed that stripes might be a common 
C~~~- ' origin of incommensurate low-energy spin fluctuations [L^^ . This view has been fueled by the 

^D I observation of stripe-like modulations in the charge sector on the surface of Bi2Sr2CaCu208+5 

^^ ' and Ca2-a;Naa;Cu02Cl2 using scanning tunneling microscopy (STM) [7J|8l[9]. It has been ar- 

gued that a picture of spatially disordered (i.e. slowly fluctuating or disorder-pinned) stripes is 
broadly consistent with both neutron and STM data [TUlIll) . In this scenario, modulations in 
the charge sector are driving incommensurabilities in the spin sector. This view is compatible 
, with the phase diagrams of various 214 cuprates, where charge order is established at a higher 

jrt ' temperature than spin order, i.e., charge stripes can exist without ordered magnetism over a 

range of temperatures. 

Recent experimental data on de- twinned YBCO, where both neutron scattering [4 and 
thermoelectric transport ;5j data show distinct in-plane anisotropies below a doping-dependent 
onset temperature, have suggested to consider electron-nematic order, i.e., rotation symmetry 
breaking with translation symmetry preserved, as a separate player. In fact, simple model 
calculations for the anisotropy of the spin excitations [T^] and that of the Nernst signal [13] 
based on uniform distortions describing nematic order were found to describe salient aspects 
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of the dataj^ This naturally prompts the question whether nematic or stripe order should 
be considered as the "primary" order parameter. This issue is of importance, as it pertains 
to proposals which link the origin of the pseudogap in underdoped cuprates to a symmetry- 
breaking order competing with superconductivity. 

In this paper, we shall attempt to answer the question whether the neutron data of Ref. [1] 
are consistent with a scenario of spatially disordered charge stripes as advocated in Ref. [TU] . 
To this end, we study the coupled order-parameter field theory of Ref. |T0] in a regime of 
a small spin gap under the influence of lattice anisotropics. Using plausible assumptions for 
the temperature dependence of the collective charge degrees of freedom, we determine the 
anisotropy of the low-energy spin fluctuations and compare them to the data of Ref. [4] . Doing 
so, we find the answer to the above question to be a tentative "yes" ; hence, we provide here 
a mechanism for the observed temperature-dependent incommensurability [J alternative to 
the one of Ref. |12j . Finally, we propose further experiments to discriminate between available 
theoretical scenarios. 

The body of the paper is organized as follows: To set the stage, we shall begin with general 
remarks about order parameters, symmetries, and phase diagrams in Sec. [H We then turn 
to the spin dynamics and describe in Sec. [3] the coupled order-parameter theory of Ref. [lOj . 
with appropriate modifications. Sec. S] presents the numerical results for the dynamic spin 
susceptibility, together with the analysis of the low-energy incommensurability. A discussion of 
experimentally relevant issues and an outlook will close the paper. 

2 Stripe and nematic order 

In the context of uni-directional stripe order in quasi-2d systems, order parameters for spin 
and charge density waves as well as for rotational symmetry breaking can be considered. In 
the following we shall discuss a continuum description which is appropriate for slowly varying 
order-parameter fields. For reasons explained below, the calculation in Sec. [3] will instead use a 
lattice formulation. 

The charge density wave (CDW) requires a pair of complex scalar fields (j)cx, (f'cy for the two 
CDW directions with wavevector Qcj; and Qcj^O Then, the charge density follows 

(p(R, r)} = pavg + Re [e^Q-^</),(R, r)] . (1) 

Similarly, there is a pair of complex vector fields (fisax , 4'say describing spin density waves (SDW) 
with wavevectors Qsa; and Qs^. For cuprates at dopings above 5%, order has been found at 
Qsx = 27r(0.5± 1/M,0.5), Q^^ = 27r(0.5,0.5 ± 1/M) and Qc:^ = (27r/7V,0), Qc^ = (0,27r/7V), 
where M and N are the real-space periodicities which follow M = 2N to a good accuracy [TJ 
[21I3]- Finally, there is an Ising scalar 0„ for / = 2 spin-symmetric electron- nematic order at 
wavevector Q = 0. The dynamics of any of these order parameters should be described by 
an appropriate 0"* (or Landau-Ginzburg- Wilson) action. Other order phenomena, e.g., uniform 
and modulated superconducting pairing, shall not be of interest here. 

Biquadratic local couplings between all order parameters are generically symmetry al- 
lowed. For instance, a term z;|(^cx(i", t)\'^\4'cv{i^, t)]'^, which couples horizontal and vertical charge 
stripes, decides about repulsion or attraction between 4>cx,y, i-e., w > will lead to uni- 
directional (stripe) order where v < results in bi-directional (checkerboard) order. 

More interesting are couplings involving one order parameter linearly. Those terms are 
strongly constrained by symmetry and momentum conservation. A nematic order parameter 
(pn in a tetragonal environment couples to CDW and SDW according to 

Hl(t>n{\4'cx\'^ - \4>cy\'^) + '«2</>n ( |0sax P " \4'say\'^), (2) 



^ Note that the YBCO crystal structure is orthorhombic due to the presence of CuO chains, such 
that the in-plane rotation symmetry is broken from the outset. Hence, there is no symmetry-breaking 
electronic phase transition, but the common assumption is that the structural anisotropics are enhanced 
by correlation effects at low temperatures, signaling the tendency toward electron-nematic order. 

^ This continuum description cannot capture modulations with wavevectors significantly different 
from Qcx and Qcy. 
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Fig. 1. Schematic phase diagrams, illustrating possible transition scenarios into symmetry-broken 
states. The vertical axis is temperature, whereas the horizontal axis - where x may represent doping - 
tunes the interplay between superconductivity and spin/charge/nematic order, which are assumed to be 
competitive. The solid (dashed) lines are transitions into lattice-symmetry-breaking (superconducting) 
states. Case a) corresponds to a weak-coupling scenario, with a single stripe ordering transition, while 
b) may be realized at strong coupling, with distinct transitions for nematic, charge, and spin order [15] . 
Case c) has a nematic, but no CDW phase: in this scenario, the nematic order may be stripe-free and 
eventually trigger a transition to an incommensurate SDW phase. In all cases, quenched disorder will 
smear out both CDW and nematic transitions and likely turn the SDW into a cluster spin glass [18| . 
Further, a structural anisotropy will smear out the nematic transition - this applies to the LTT phase 
of 214 cuprates and to YBCO. Qualitatively, case b) might be realized in some 214 cuprates (with the 
nematic transition replaced by the LTO-LTT transition), whereas cases b) or c) might apply to YBCO. 



note that |0ca;P etc. carry vanishing lattice momentum. A CDW order parameter couples to a 
spin density wave (/)« according to 

'^sl^c^L + c.c.) 



(3) 

which is allowed only if the ordering wavevectors obey Qc = 2Qs. In addition, symmetry- allowed 
higher-order couplings between gradients of order parameters will occur. 

The terms (J2l3p imply that a uni-directional density wave induces nematic order via k\^i 
and that a coUinear SDW with Qg induces a CDW with 2Qs via K3. As a consequence, a 
transition into a stripe-ordered state may occur as a direct transition from a disordered into 
a CDW+SDW state, or via intermediate nematic and CDW phases, see Fig. [T] The phase 
diagram of the corresponding Landau theory has been worked out in detail by Zachar et al. [1A\ . 
Whether the intermediate phases are realized depends on microscopic details. In particular, an 
electron-nematic phase can come in (at least) two flavors, namely (i) a purely uniform (Q = 0) 
distortion, as arises from a Pomeranchuk instability of the Fermi surface (dubbed "stripe-free" 
nematic in the following), or (ii) a precursor to a uni-directional stripe phase, in which case 
it will display substantial and slow charge fluctuations at non-zero wavevectors (hence dubbed 
" stripe- driv en") . Although stripe-free and stripe-driven nematics do not differ in symmetries, 
the underlying microscopic mechanisms are clearly distinct. 

While weak-coupling theories typically give a direct transition into a CDW-I-SDW state, 
Fig. [T^, distinct scenarios are viable at strong coupling [15j . For instance, the tendency toward 
bond order could be dominant [16j . which induces both translation and rotation symmetry 
breaking in the charge sector. Incommensurabilities in the spin sector would then follow from 
periodic modulations of the charge, as in Fig. [T}d. Alternatively, stripe-free electron-nematic 
order might be the primary phenomenon, which then drives ordering phenomena in both the 
charge and spin sectors. Fig. [1}:. 

For 214 cuprates with LTT lattice structure, the in-plane rotation symmetry is broken in 
each Cu02 layer, such that there can be no transition into a purely nematic phase. Experimen- 
tally, distinct transition are seen in the charge and spin sectors in La2-a:Baa;Cu04 as well as in 
La2-a;Sra;Cu04 co-doped with Nd or Eu. These transitions can be associated with the onset of 
CDW and SDW order, respectively. Remarkably, these transitions apparently continue to ex- 
ist in Lai5/8Bai/8Cu04 under high pressure where the material remains in the so-called HTT 
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phase, without structurally broken in-plane rotation symmetry, down to lowest temperatures 
[17]. In contrast, in YBCO, no signatures of (static) charge order have been reported. Static spin 
order is restricted to underdoped samples and very low temperatures: in YBCO-6.45, quasi- 
static incommensurate order in neutron scattering sets in below 30 K, with a correlation length 
of about 20 A only, and /xSR measurements detect static moments only at 1.5 K, suggesting 
that the order is weak and glassy [3]. However, signatures of rotation symmetry breaking are 
prominent in both the spin excitations W and in the Nernst effect 5 . Although YBCO also 
displays a structural anisotropy, smearing out a transition into an electron-nematic state, the 
experimentally observed anisotropics display a relatively well-defined doping-dependent onset 
temperature which has been proposed to follow the pseudogap temperature. 

At present, it is unclear whether the nematic order in YBCO should be viewed as stripe- 
free or stripe-driven. In the literature, different points of view have been put forward. On the 
one hand, simple stripe-free models of electron-nematic order have been invoked to explain 
the anisotropic incommensurate spin excitation spectrum in YBCO as detected in Refs. [H 
\W\ . Here, both RPA on top of an anisotropic band structure j20llMl l22 , 23 and a Ginzburg- 
Landau theory for spin excitations with phenomenological anisotropy ^Tl^ have been employed. 
A model of electrons moving in an anisotropic band structure has been found [13] to describe 
some features of the Nernst effect measured in field up to 15 T in YBCO [Sj- The Nernst analysis 
also suggests that stripe order sets in at lower temperatures [Mll25| . On the other hand, the 
notion of apparently "universal" spin excitations at intermediate energies in various cuprates 
[6] has triggered stripe-based explanations. In particular, fluctuating or disorder-pinned stripes 
have been shown to induce a hour-glass-like spin excitation spectrum which is compatible with 
most of the YBCO neutron scattering data. In this scenario, anisotropics in x" arise due to 
the small (chain-induced) bare electronic anisotropic which is strongly enhanced due to the 
proximity to a uni-directionally ordered state. 

3 Lattice order parameter theory for disordered stripes 

The remainder of the paper is concerned with the question of whether a scenario of disordered 
stripes can be consistent with the neutron data obtained on YBCO-6.45 [J. In this section, 
we summarize the coupled order-parameter theory of Refs. 10,26J and discuss the parameter 
regime of potential relevance for the experiment. Numerical results follow in the next section. 

The approach of Refs. [TIJ[[2S] builds on a few main assumptions: (i) The collective physics 
in both the spin and charge sectors can be captured by order-parameter fields, without ref- 
erence to single-particle excitations, (ii) The "bare" spin sector is dominated by fluctuations 
near wavevector (tTjTt), and any incommensurability in the spin sector is therefore driven by 
its coupling to the charge sector, (iii) In the presence of incommensurabilities, the spin order 
remains coUinear. Item (ii) implies that neither a continuum formulation nor a perturbative 
treatment of the spin-charge coupling are justified, because the spin ordering wavevector can 
change significantly (and in a non-perturbative fashion) due to a coupling to the charge modes. 
Items (i) and (ii) are borrowed from the successful description of magnetic excitations of static 
stripes, as measured by Tranquada et al. [57] in stripe-ordered Lai5/8Bai/8Cu04, using mod- 
els of coupled spin ladders [55[f^[501l5T] . Such models assume static charge order, which in 
turn imposes a strong modulation on both spin density and magnetic couplings; without this 
modulation the spin sector is described by a nearest-neighbor Heisenberg model. 

For the lattice order-parameter description, we start from a (^^ theory for the spin sector: 

S^ - JdTj2 [idrV^jf + ^^?] + J2''^i^j-^ff + Si (4) 

J ijj') 

with 54 being the quartic self- interaction term. The real field ipj = fi^j) is defined on the 
sites j of the square lattice of Cu atoms. It carries lattice momentum Q — {tt,tt), such that 
(fj and the Cu spins Sj are related through Sj oc e**^'"^(/9j, and ipj = const, describes the 
familiar two-sublattice antiferromagnet. In the absence of static order, the quartic term 54 may 
be neglected, which amounts to a Gaussian approximation to the spin fluctuations. 
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For the charge sector we introduce complex order parameters 'ipx,y for stripe modulations 
with wavevectors Qcx and Qcj, - below the focus will be on period-4 stripes. Then, the real 
field Qx{j^) = Ref/'a:(r)e**^"^'' (similarly for Qy) measures the modulation of site and bond 
charge densities. We choose signs such that Sp{rj) = Qx + Qy is the deviation of the local hole 
density from its spatial average. The allowed couplings between spin and charge are dictated by 
spin and time reversal symmetries, with the lowest order ones being of the form XQ(p^ . Guided 
by the lattice models [2511^ we choose [2S] 



S^^ = / dr ^^XiQ xir j)(p'^j + \2Qx{rj+x/2)'P J 



V]+x 



+ X3Qx{''^j)fj~x'f3+x + MQx{'rj+y/2)Vj^3+y 

+ X -^ y 



(5) 



Here, Ai > implements the correlation between the on-site charge and spin densities, i.e., a 
large hole density reduces the local amplitude of the spin fluctuations. A2-4 ensure that the 
effective first- and second-neighbor exchange constants modulate along with the charge/bond 
order. In particular, the antiphase domain wall properties of the stripes |32[|33) are refiected in 
the positive sign of A2,3, e.g., for large \2Qx the nearest-neighbor magnetic coupling perpendic- 
ular to the stripe direction in S^p -\- S^^ becomes ferromagnetic (instead of antiferromagnetic) . 

Perfect static vertical (horizontal) stripe order corresponds to ^x — const, ipy — (ipy = 
const, ipx = 0)- Then, the action S^p + S^^ is a theory for magnetic modes in a background of 
static charge order. For sufficiently large A couplings, the minimum energy of the (/? fluctuations 
will be shifted away from (tt, tt) to the incommensurate wavevector Qs dictated by the charge 
order (which fulfills 2Qs = Qc). A Gaussian approximation to such a theory 26| reproduces 
the results obtained from microscopic models of coupled two-leg ladders ^28.29^ which are in 
agreement with experimental data on Lai5/gBai/gCu04 [77] . 

Deviations from perfect charge order can be described by a f/'^-type theory of the form 



S,h = / dT^Udri'ixf' + \drll^iyf' -t- Sx\'4'ix\'^ + SylV'jyP 



+ c\^\'4'%x - i^i+x.x? + clxl^ix - i'i+y,x? + c\y\lpiy - i>t+x,yf + cly\lpiy - 1pi+x.y\^ 

+ u.rPt + u^i^t + v\i^.x?\^,y? + w {^l+^*^+^l + ^*^)\ (6) 

with %j:ix = ^pxiY^) and ipf = \ipix\'^ + \''Piy\'^- This theory features a variety of non-linear terms, 
which have been included to model the non-trivial physics of the underlying strongly-correlated 
system on the lattice scale (recall that the objective is not to describe universal critical phe- 
nomena). Here, a combination of ui < and U2> suppresses amplitude fluctuations of f/; - 
this is motivated by the STM data of Ref. "Q". For cix = ciy, C2x — C2y, Sx — Sy, and v — w — 0, 
the action has 0(4) symmetry, which, however, is broken by various lattice effects: The quartic 
ti|'!/'a;P|''/'y P term regulates the repulsion or attraction between horizontal and vertical stripes; 
we employ ■i;>0 leading to stripe-like order (whereas ti<0 results in checkerboard structures). 
The phase-sensitive w term provides commensurate pinning and selects bond-centered (instead 
of site-centered) stripes pfTO] for w>0. 

The treatment of 5^-1-5^-1- 5,^^ (01-0 requires further approximations. Ref. J^ proposed 
an adiabatic approximation, based on the assumption that the charge modes are slow compared 
to the spin modes. As in the familiar Born-Oppenheimer approximation, the spin sector may 
then be treated for a fixed configuration of charges, and a proper average over the charge 
configurations has to be taken. For simplicity, one may further neglect the feedback of the spin 
sector on the charge sector. Practically, this then amounts to the following procedure: A classical 



^ Linking the notations of Sec. [2] and Ref. [TD], (j)c = ^ and (j>s = (p. However, -tp and (p are lattice 
instead of continuum fields. 
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Monte-Carlo simulation of the action S.^i, (without the dr terms, at a temperature T = 1) is used 
to generate individual configurations of the charge order parameter field ipx,y (on a finite lattice 
of up to 64^ sites), from which one calculates the Qx.y For each configuration of Qx,y, S^ -\-S^^ 
defines the action for the spin sector. Neglecting the quartic term ^4 in Eq. (|1]), S^p + S^^ can be 
diagonalized numerically. The resulting spin susceptibility x" is finally averaged over typically 
20 Monte-Carlo configurations. While this procedure misses all inelastic scattering processes, it 
properly captures the influence of short-range quasi-static stripe order on the magnetic modes. 
Clearly, it is appropriate both for slowly fluctuating and for static but spatially disordered 
stripes - the latter is expected from impurity pinning which acts as a random field on the 
CDW order parameter. 

Although the theory iS^ + iS^ + iS^^ has a number of input parameters, most can be fixed 
with an eye towards experiment, c in iS^ sets the velocity of spin excitations and is known 
experimentally; s is then chosen as to reproduce the energy of the spin excitations at (tt, tt) 
(the "spin resonance"). The overall scale of the parameters in <S^ sets the typical amplitude 
of the charge fluctuations IV'ltyp- The product Aji/jltyp enters S^p^p and is chosen such that the 
spin excitations become gapless at the incommensurate wavevector Qs for perfect stripes. The 
interplay of the remaining parameters ci^2;x,yi Sx,y, ui^2, v, and w determines the character of 
the stripe fluctuations. Here, reasonable combinations can be found be comparing the simulation 
result for Sp — Qx + Qy with charge configurations deduced from STM data [9J. The latter 
indicate that amplitude fluctuations are suppressed as compared to phase fluctuations, that 
bond-centered pinning is strong, and that domain walls are rather sharp. Finally, keeping all 
other parameters fixed, a variation of Sx^y will drive the charge sector through its order-disorder 
transition. We note that the theory 5^ ^ can be expected to have a nematic phase at large 
positive v; however, we have not mapped out the full phase diagram of S^. 

Interestingly, Ref. [TU] first determined plausible sets of parameters which were able to 
reproduce an hour-glass spin fluctuation spectrum within the adiabatic approximation to S^ + 
S^ + Sip^ in a regime where the charge correlation length was between 10 and 50 lattice 
constants. Later, the action 5^, coupled to the single-particle sector of a d-wave superconductor, 
was used to model electronic spectra of disordered-pinned stripes. Here, the same parameters 
for <S^ were used and good agreement with the STM data of Ref. [5] was obtained. 

Spatial anisotropies arising from structural distortions of the Cu02 layers enter the theory 
at various places. Most important are the parameters of iS^, in particular if the charge sector 
is close to spontaneous rotational symmetry breaking and hence sensitive to small anisotropies. 
Below, we shall employ anisotropic mass terms, Sx 7^ Sy, but anisotropic gradient terms are 
feasible as well. 

We close this section with a remark on Ref. |12) . which proposed a mechanism for incom- 
mensurate spin excitations induced by uniform nematic order. In fact, this mechanism bears 
similarity to the one implemented in S^ + S^^ used here and in Ref. [10] : Using the notation of 
Ref. [12] , the gradient term in Eq. (4) of Ref. fT^ has a negative prefactor in one direction for 
C2N > Cq - the same happens in S^ +Sip^ for large positive A2. Microscopically, this translates 
into effectively ferromagnetic (instead of antiferromagnetic) coupling, and higher-order gradi- 
ents are required to render the theory well-defined. However, from Ref. [12j it remains unclear 
why uniform nematic order should induce ferromagnetic effective couplings. Also, for the ex- 
perimentally relevant ordering wavevectors, the gradient expansion of Ref. [12] is questionable, 
and a lattice theory appears more appropriate. 

4 Spin excitations in the presence of lattice anisotropy 

We now present numerical results for the magnetic excitation spectrum of disordered stripes 
in the presence of lattice anisotropies. Although initial results were reported in Ref. [10], those 
were in a regime of a sizeable spin gap, and no attempt was made to model the temperature 
evolution. In contrast, here we will choose parameters with an eye towards the experimental 
data on YBCO-6.45 of Ref. [4] where the spin gap is small [34]. Such variations of the spin gap 
are illustrated in Fig. [5] where we show results for x"(q, w) with varying mass parameter s in 
S^, here without in-plane anisotropy. 
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Fig. 2. Neutron scattering intensity, x"(q,'^) for momentum {qx,n), from Monte-Carlo simulations 
of the coupled order parameter theory, Eq. (|2^, for varying spin-gap parameter s in S^p |34II35| . 
The energy has been normalized to that of the "resonance" peak (which is roughly 35 meV for the 
employed parameters [36]) and an artificial energy broadening of -T = 0.05 has been used. Parameters 
are s^ ~ Sy = —3.4 and Mi = —1.15, placing S^ somewhat inside the charge-ordered phase, and 
a) s = 105 meV, b) s = 88meV, c) s = 70meV [36) . There is no in-plane anisotropy in the input 
parameters, and the spectra have been symmetrized w.r.t. q^ ^ Qy The hour-glass shape is clearly 
visible. 



To quantify the incommensurability of the spin excitations, we shall analyze the low-energy 
cuts through x"(<li ^) in the presence of an in-plane anisotropy in a manner similar to Ref. [1] in 
order to extract an apparent incommensurability of the spin fluctuations. To assess the changes 
in the spin-fluctuation spectrum with temperature, we shall vary the input parameters for S^ 
such that both amplitude and correlation length of stripe fluctuations change. 

The following assumptions on parameters will be made: (a) The spatial period of the incip- 
ient charge order (and thus Qcx, Qcy) will be held fixed - this is motivated by the essentially 
temperature independent stripe modulation period in 214 cuprates. For reasons of compu- 
tational simplicity and momentum resolution, we shall choose a charge period of 4 lattice 
constants. (Note that period 8 might be more appropriate for YBCO-6.45; as a result, our 
incommensurability of the spin excitations will be roughly twice as large as in experiment.) (b) 
We shall vary the strength of the charge order by varying Sx,y and ui - these parameters de- 
termine the correlation length and mean amplitude of i'x,y Qualitatively, these changes mimic 
the influence of temperature on the charge sector. (Recall that the simulation of 5^ is classical.) 
Other parameters will be held fixed, (c) The lattice anisotropy of YBCO will be implemented 
by anisotropic mass terms in S^, Sx j^ Sy. Near the ordering transition in the charge sector, 
this will induce strongly anisotropic charge correlations, which in turn lead to anisotropics in 
the spin sector. (Alternatively, anisotropic velocities in 5^ can be employed.) In principle, a 
small anisotropy should also be present in the bare spin action S^, but this will be neglected 
as the spin anisotropy is dominated by the influence of the charge sector. 

A central question is whether incipient charge order with a fixed spatial period can induce 
incommensurate spin excitations with a temperature-dependent apparent incommensurability. 
Ref. [H] suggested that this is not the case (based on a mean- field argument only) , which then 
might rule out incipient stripe order as explanation for the observation in Ref. |3|. 

A sequence of results for x"(qjw) with varying strength of the charge order is shown in 
Fig. [3| For the parameters in panels b) and c), the correlation length in the charge sector is 10 
or smaller, and the coupling to the spin sector produces a broad spin excitation spectrum at 
low energies, with a "vertical dispersion" , very similar to observations in underdoped YBCO 
[lirrS]. Note that moving closer to or inside the charge-ordered phase produces an hour-glass 
spectrum in x"(q, w) (see Figs. [Hand [3^ as well as Figs. 2,3 of Ref. [TU]). 

The data in Fig. [3] have been calculated in the presence of a small in-plane anisotropy. 
Consequently, x"(q, w) at constant uj below the "resonance" energy is anisotropic along the x 
and y directions, middle panel of Fig. [31 The anisotropy at higher energies is less significant 
(not shown). Following the experimental analysis of Ref. [3], we can fit the broad intensity 
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Fig. 3. Results for x"(cii'^) |34II35] in the presence of an in- plane anisotropy for varying stripe 
correlation length. Top: Color-coded intensity as function of energy and momentum ((7a:,7r) along the 
X axis. Middle: Constant-energy scans aX u = 0.4. Bottom: Intensities along (qx,ii) and {Tr,qy) at 
oj = 0.4, showing the low-energy anisotropy. The evolution from panels a to c mimics the increase of 
temperature: Here, both Sx,y and ui in the action S^ for the charge order parameter are varied, such 
that the stripe correlation length ^ decreases from a) to c). The anisotropy is encoded in the mass 
difference for horizontal and vertical stripes, (sx — Sy), and is held fixed, as is the spatial periodicity, 
A'^ = 4, of the underlying charge stripes. Parameters are Sy — s^ = 0.01, a) Sx ~ —2.8, ui = —1.13, 
5 ^ 12, b) Sx = -2.6, ui = -1.11, C ~ 9, c) s^ = -2.0, m = -1.05, ^ ^ 4 [^. 



profile along the x direction, bottom panel of Fig. [3l to a sum of two Gaussians, centered at 
7r(l ± es, 1). The apparent incommensurability e^ resulting from such a fit is plotted in Fig. |3ti 
as a function of the mass of the charge order parameter. A number of features can be noted: 
Cs is finite both inside and outside the charge-ordered phase. It effectively vanishes once the 
stripe correlation length becomes shorter than 5 lattice spacings (this value depends on various 
parameters in <S^ -I- 5^^), and it approaches its "saturation" value es = l/N dictated by the 
charge order only deep in the ordered phase. Hence, even long-range ordered stripes do not 
necessarily result in an apparent incommensurability of es = 1/A^. 

Fig- Eh- niakes clear that the apparent incommensurability e^ obtained from fitting the profile 
of finite-energy spin fluctuations cannot be directly used to deduce the spatial period of possibly 
underlying stripe order. In particular, a temperature dependence of Cg as observed in Ref. [4] 
does not imply a temperature dependence of the underlying stripe period, but can simply be a 
result of smearing due to thermal and quantum fluctuations. 
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Fig. 4. a) Apparent incommensurability Es of the low-energy intensity profile in x"(q,'^) as function 
of the order-parameter mass s^; ui is varied along with s^ according to ui = Sx/W — 0.85. The two 
curves correspond to Sy — s^ — 0.01 (solid) and 0.02 (dashed); with the charge-ordering transitions 
located at Sx values of roughly 2.9 and 2.8, respectively, e^ is obtained from fitting x" at a; = 0.4 along 
{qx, tt) to a sum of two Gaussians centered at 7r(l ± e^, 1). (ej, is set to zero when the fit with a single 
Gaussian becomes visually indistinguishable from the fit with two Gaussians.) For perfect stripes, the 
incommensurability is ts = 1/4. b) Momentum-integrated structure factor S{iu) as function of energy 
for different stripe correlation lengths. The inset shows the structure factor as obtained from the 
results for x"(ci, i^), the main panel shows the same data with a "background" subtracted - here the 
reference curve is chosen as the one with the shortest correlation length in the inset. The background 
subtraction mimics the corresponding experimental procedure, where the signal at high temperatures 
is typically smeared out due to a combination of thermal fluctuations and fermionic damping and a 
linear background is subtracted. (The noise in all curves is due to Monte Carlo sampling.) Parameters 
are Sy — Sx = 0.01, Sx = —2.8, wi = —1.13 (solid), Sx — —2.4, ui = —1.09 (dashed), Sx = —2.0, 
Ml = —1.05 (dash-dot), Sx ~ —1.0, ui — —0.95 (short dash) [36) . 



Finally, we show in Fig. |4}d the momentum-integrated dynamic structure factor S{io) ob- 
tained from the numerical simulations |35) . One observes a clear transfer of spectral weight 
from high to low energies upon with increasing charge correlation length, not unlike in the 
experimental data [1]. Note that, in the present data, S{uj) never displays a maximum at low 
energies (inset of Fig. HId); instead, a plateau emerges near the resonance energy which turns 
into a peak for larger stripe correlation length. A low-energy peak in ^(a;) can only be obtained 
after a "background" subtraction (main panel of Fig.HjD), here a "high-temperature" curve has 
been subtracted. In the experiment, background subtraction is done as well, with the difference 
that the high-temperature signal has much less structure due to thermal spin fluctuations and 
fermionic damping, both of which are absent from our simulation. 



5 Conclusions 



This paper has employed a scenario of slowly fluctuating or disorder-pinned stripes to model the 
anisotropic spin excitation spectrum of underdoped YBCO. We have found that such a theory, 
using plausible input parameters, is able to reproduce salient features of the experiment of 
Ref. [4], including the temperature-dependent apparent incommensurability in the low-energy 
spin excitations, in contrast to what has been suggested in the literature. 

Therefore, both stripe-driven and stripe-free nematic orders appear viable candidates to 
explain the low-temperature anisotropics in YBCO. It remains to ask which experiments could 
distinguish the two scenarios. Clearly, a direct observation of finite-q dynamical charge fluc- 
tuations would support the stripe scenario, but a suitable experimental probe does not exist. 
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In the presence of quenched disorder, stripes will be pinned and should be observable in real 
space, e.g., by STM. Unfortunately, clean surfaces as required for STM are hard to obtain for 
this cuprate family. Resonant X-ray scattering can only detect static charge order if the corre- 
lation length is sufficiently large, which is probably not the case in YBCO [5]. A more indirect 
probe is via the Fermi surface geometry. Whereas a stripe-free nematic leads to a distorted 
single-sheet Fermi surface (at least in the simplest tight-binding models), a stripe phase pro- 
duces, depending on the periodicity, a variety of smaller pockets [37,38 . Thus, photoemission 
experiments with sufficient energy resolution might eventually be able to distinguish the two 
scenarios. (Note that for fluctuating stripes those pockets will not be present at lowest ener- 
gies; this is also of relevance for the interpretation of Nernst data, where it has been argued 
[24) that static stripes are incompatible with the large Nernst anisotropy seen at intermediate 
temperatures [S].) Present-day quantum oscillation data point toward the existence of pockets, 
albeit in a large magnetic field [39] . 

I thank R. K. Kaul, S. Sachdev, T. Ulbricht, T. Vojta, and A. Wollny for collaborations and fruitful dis- 
cussions on this subject, E.-A. Kim, M. Lawler, and K. Sun for correspondence regarding Ref. 12 , and 
furthermore E. Fradkin, M. Greiter, V. Hinkov, B. Keimer, S. Kivelson, J. Tranquada, and G. Uhrig for 
discussions and correspondence. This research was supported by the Deutsche Forschungsgemeinschaft 
through the Research Unit FOR 538. 
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